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We study a dynamical Ising model of agents' opinions (buy or sell) with coupling coefficients 
reassessed continuously in time according to how past external news (magnetic field) have explained 
realized market returns. By combining herding, the impact of external news and private information, 
we test within the same model the hypothesis that agents are rational versus irrational. We find 
that the stylized facts of financial markets are reproduced only when agents are over-confident and 
mis-attribute the success of news to predict return to herding effects, thereby providing positive 
feedbacks leading to the model functioning close to the critical point. 
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PACS numbers: 64.60.Ht; 89.65.Gh; 87.23.Ge 



Social systems offer a fascinating field for the applica- 
tion of recent concepts and methods developed in Physics 
to tackle complex ./V-body systems with nonlinear feed- 
backs, and many competing states. A long tradition 
started with the application of Ising models and its exten- 
sions to social interactions and organization 0, 0, 0, • 
A large set of economic models can be mapped to various 
versions of the Ising model to account for social influence 
in individual decisions (see Q and references therein). 
Other recent works usingthe Ising model include models 
of bubbles and crashes Q , a version with stochastic 
coupling coefficients which leads to volatility clustering 
and a power law distribution of returns at a single fixed 
time scale 8] , and models of opinion polarization [9|, llOj . 
The dynamical updating rules of the Ising model can be 
shown to describe the formation of decisions of boundedly 
rational agents or to result from optimizing agents 
whose utilities incorporate a social component [5;]. The 
Ising model is one of the simplest models describing the 
competition between the ordering force of imitation or 
contagion and the disordering impact of private informa- 
tion or idiosyncratic noise, which leads already to the 
crucial concept of spontaneously symmetry breaking and 
phase transitions |l2( . 

However, human beings are not spins, they can learn, 
that is, adapt the nature and strength of their interac- 
tions with others, based on past experience. In the lan- 
guage of the Ising model, this amounts to generalize to 
time-dependent coupling coefficients which reflect past 
experience. Here, we study a generalized Ising model of 
interacting agents buying and selling a single financial 
asset who base their decisions on a combination of mu- 
tual influences or imitation, external news and idiosyn- 
cratic judgements. The novel ingredient is that they up- 
date their willingness to extract information from the 
other agents' behavior based on their assessment of how 



past news have explained market returns. Agents up- 
date their propensity to herding according to how the 
news have been successful in predicting returns. We dis- 
tinguish between two possible updating rules: rational 
and irrational. In the rational version, agents decrease 
their propensity to imitate if news have been good pre- 
dictors of returns in the recent past. In the irrational 
version, agents mis-attribute the recent predictive power 
of news to their collective action, leading to positive self- 
reinforcement of imitation. We show that the model can 
reproduce the major empirical stylized facts of financial 
stock markets only when the updating of the strength 
of imitation is irrational, providing a direct test and the 
evidence for the importance of misjudgement of agents bi- 
ased toward herding. This model also offers a dynamical 
derivation of the multifractal properties of the structure 
functions of the absolute values of returns and their con- 
sequences in the characteristic power law relaxations of 
the volatility after bursts of endogenous versus exogenous 
origins. 

Consider n 2 agents interacting within anxn 2D square 
lattice network J\f (we have verified that the properties 
described below are not sensitive to n in the range 20 — 
100 that we tested). At each time step t, agent i places 
a buy (si(t) = +1) or sell (sj(t) = —1) order according 
to the following process 



Si(t) = sign 



E 



i^(t)E[ Si ](f) + *i(t)G(t) + ei(t) 



(1) 

where E[sj](i) is the expectation formed by agent i on 
what will be the decision of agent j at the same time t. 
An agent i imitates only her friends, that is Kij = if 
i = j or j is not a friend of (connected to) i. Expres- 
sion Q embodies three contributions: (i) Imitation in 
which is the relative propensity of the trader i to 
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be contaminated by the sentiment of her friend j; (ii) 
The impact G(t) of external news (positive resp. nega- 
tive for favorable resp. unfavorable news) and CTj(£) is the 
relative sensitivity of agent's sentiment to the news, uni- 
formly distributed in the interval (0, <7 max ) and frozen to 
represent the heterogeneity of the agents; (iii) Idiosyn- 
cratic judgement ej(i) associated with private informa- 
tion, assumed to be normally distributed around zero 
with an agent-dependent standard deviation s e< i equal 
to the sum of a common constant CV and of a uniform 
random variable in the interval [0,0.1] again to capture 
the heterogeneity of agents. We have tested several im- 
plementation of the formation of expectations E[sj](t) in 
JU, such as backward looking (Ej[sj](t) = Sj(t—1) for all 
i and j's) or information cascades along specific chains 
within the network which give similar results. 

We account for the adaptive nature of agents and their 
learning abilities by updating the coefficient of influence 
of agent j on agent i according to the following rule: 

Kij{t) =Ki(t) =b i + aK i (t-l) + (3r(t-l)G(t-l) . (2) 

The idiosyncratic imitation tendency hi of agent i is uni- 
formly distributed in (0,6 max ) and frozen. The coeffi- 
cient a quantifies the persistence of past influences on 
the present (the case a = has been studied in part in 
Ref. It captures the fact that social connections 

evolve slowly and exhibit significant persistence, as doc- 
umented in numerous studies |l4| . Networks of investors 
communicating their opinions and sentiment on the stock 
market are similarly persistent. The last term with (3 ^ 
quantifies how agent i updates her propensity for imita- 
tion based on the role of the exogenous news G(t) in de- 
termining the sign and amplitude of the observed return 
in the preceding time period. Note that, by construction, 
the model is non- variational as Kji(t) ^ Kij(~t) in gen- 
eral if the biS arc different, but this is not crucial for our 
results. Finally, the market price is updated according 
to p(t) — p(t — l)exp[r(i)] where the return r{t) obeys 

r W = ^ ie \N ; where A measures the market depth or 
liquidity. The return is thus proportional to the "mag- 
netization" or aggregated decisions of the agents. 

The sign of the coefficient j3 is crucial. For j3 < 0, agent 
i is less and less influenced by other agents, the better has 
been the success of the news in determining the direction 
and amplitude of the market return. This process is self- 
reinforcing since, as K+j decreases, the dominant term 
becomes <Ji(t)G(t), which further ensures that the news 
correctly predict the decision of agents and therefore the 
direction of the market move, thus decreasing further the 
coefficient of influence Kij . Reciprocally, agents tends to 
be more influenced by others when the news seems to in- 
correctly predict the direction of the market. The news 
being not reliable, the agents turn to other agents, be- 
lieving that others may have useful information. This 
is in agreement with standard economics which views 



the stock market as an efficient machine transforming 
all news into prices. 

For j3 > 0, the more the news predict the direction of 
the market, the more the agents imitate other agents. 
This is the "irrational" case where agents either mis- 
attribute the origin of the market moves to herding rather 
than to the impact of news, or misinterpret the exogenous 
character of news in terms of endogenous herding or in- 
fer that other agents will be following more eagerly as a 
group the direction given by the news. This may occur 
due to mutually-reinforcing optimism [l5| and overconfi- 
dence 

In our simulations, we fix A = 40 to obtain returns with 
amplitudes comparable to that of empirical observations 
and a = 0.2. Similar results are obtained for a — 0.4, 0.6, 
and 0.8. We have explored the properties of the model in 
the parameter space of b max , cr max and CV. There is no 
loss of generality in fixing \(3\ = 1 to explore the relative 
importance of the term /3r(t— l)G(t— 1), since the typical 
scale of the Ki's is set by 6 max whose amplitude is varied 
in our numerical exploration. 

For j3 = —1, it is easy to show that the attractor of 
the dynamics is characterized by negligible imitation and 
only the news and private information terms are impor- 
tant for the dynamics. Indeed, starting from large Ki(t)'s 
such that the system is above its critical value K c and 
the agents behave collectively (ferromagnetic phase), the 
news G(t) acts as a "magnetic" field which orders the 
agents' decision accordingly, leading to the news correctly 
predicting the returns. Since f3 = — 1, the coupling coef- 
ficients i^i's will be decreased by the amount \r(t)\. This 
will continue until the K^s are in majority much smaller 
than K c , at which point the dynamics becomes stable 
in the "paramagnetic" phase because the collective deci- 
sion J^. Si and therefore the market return have little or 
no relationship with the external news. Hence, the term 
(3r(t + l)G(t) takes random signs from one time step to 
the next, leading to an effective random forcing added 
to the autoregressive equation Ki{t) = hi + aKi(t — 1). 
We thus expect Gaussian distributions of returns when 
6,;/(l — a) is smaller than K c and bimodal distributions 
when bi/(l — a) > K c reflecting the slaving of the global 
opinion to the sign of the news. Our simulations, which 
have scanned 480 different models for 6 max from 0.1 to 
0.5 with spacing 0.1, tr max from 0.005 to 0.08 with spac- 
ing 0.005, and CV from 0.1 to 1.1 with spacing 0.2, 
confirm this prediction. Consider the distribution of re- 
turns r T (t) = \n[p(t)/p(t — t)] at different time scales 
t. For large idiosyncratic noise (large CV) and not too 
large & max , the distribution of returns is Gaussian for 
all time scales r. For smaller CVs and larger & max , we 
observe multimodal return distributions. In the param- 
eter space that we have explored and notwithstanding 
our best attempts, we have not been able to find a set of 
parameters leading to distributions of returns exhibiting 
a monomodal shape with fat tails for small time scales, 
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evolving slowly towards Gaussian distributions at large 
time scales, as can be observed in empirical data |l7j . 
In addition, the correlation function of returns (C T (r, r)) 
and of volatilities (C T (|r|, |r|)) have similar amplitudes 
and decay with the same characteristic time scale as a 
function of time lag. This is very different from the ob- 
served correlations of financial markets, with very short 
memory for returns and long-memory for the volatility. 

For P = 1, we obtained the following main stylized 
facts of financial stock markets: (i) distributions of re- 
turns at different time scales r (Fig. ; (ii) correlation 
function of returns and of the absolute value of the re- 
turns (Fig. EJ); (iii) scaling of the moments of increas- 
ing orders of the absolute values of the returns (testing 
multifractality) ; (iv) the existence of a hierarchy of ex- 
ponents controlling the relaxation of the volatility after 
an endogenous shock, another hallmark of multifractal- 
ity (Fig. I2J); (v) the existence of bubbles and crashes. 
We have explored 160 models with 6 max from 0.1 to 0.5 
with spacing 0.1, <r max from 0.01 to 0.08 with spacing 
0.01, and CV from 0.1 to 0.7 with spacing 0.2. For each 
model, we generate time series of length equal to 10 5 
time steps. We have found several parameter combina- 
tions which lead to realistic stylized facts, for instance, 
(6 max , cr max , CV) equal respectively to (0.3, 0.03, 0.1), 
(0.4, 0.04, 0.1), (0.4, 0.05, 0.1), (0.5, 0.06, 0.1), (0.1, 0.01, 
0.3), (0.1, 0.02, 0.3), (0.2, 0.02, 0.3), (0.2, 0.03, 0.3), (0.3, 
0.04, 0.3), (0.5, 0.05, 0.3), (0.5, 0.07, 0.3), (0.3, 0.03, 0.5), 



and (0.5, 0.05 

(^max — 0.3, fJmax 



0.5). The results presented here are for 
0.03, CV = 0.1) which is typical. 




FIG. 1: (Color online) Probability distribution density (in 
logarithmic scales) of log-returns at different time scales r of 
a price time series obtained with A = 40, a — 0.2, 6 ma x = 
0.3, <r max = 0.03, and CV = 0.1. The log-returns r> are 
normalized by their corresponding standard deviations o>. 
The pdf curves are translated vertically for clarity. The thick 
dashed line is the Gaussian pdf. The other dashed lines have 
been obtained by r-fold convolutions of the pdf of the one- 
time step return ri(t) = \np(t)/p(t — 1). 

Figure 2] shows the evolution of the pdf 's of returns 
from stretched exponential or power laws at short time 



scales that cross over smoothly to a Gaussian law at the 
largest shown time scale, in excellent agreement with em- 
pirical facts 01- Note the difference between the con- 
tinuous and dashed lines for r ~ 4, 16, and 64, which 
expresses the existence of significant dependence in the 
time series of returns. Such behavior is very similar to 
what is observed in real data. 
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FIG. 2: (Color online) Impact of a on the auto-correlation of 
the absolute values of the returns. 

The temporal correlation of the log returns r\ as a 
function of the time lag I exhibits a very short correla- 
tion time, of duration smaller than one time step (not 
shown). In contrast, the temporal correlation of the ab- 
solute value of log returns r\ ("volatility"), taken as a 
proxy for the volatility, exhibits long-range dependence 
up to approximately 1000 time steps. The linear-log re- 
lationship suggested by the plots in Fig. [21 are predicted 
by the multifractal random walk (MRW) [18j , which pro- 
vides an excellent model of many properties of financial 
time series 0] . The MRW depends only on three param- 
eters: the multifractal parameter A 2 w 0.02 — 0.04, the 
integral time scale T « 1 year and the standard devia- 
tion of returns. Comparing the dependence properties of 
the returns and of the volatility suggests that one trad- 
ing day corresponds roughly to 5 time steps of the model. 
This correspondence translates into a integral time scale 
T of about 200 days, which is compatible with empiri- 
cal estimates for the MRW |2(|. The MRW also predicts 
(and this is well-verified by empirical data) that the au- 
tocorrelation functions of |r T (i)| for different r should 
superimpose for time lags larger than their respective r 
[18j. This prediction is also approximately observed in 
our model (not shown). 

Another important stylized facts is the multifractal 
structure of the absolute values of log-returns |2jJ, l2l| . 
We verify the existence of a strong multifractality in our 
time series (not shown) expressed by the scaling of the 
structure function M q {r) = (\r T \ q ) ~ r^ q , with exponents 
£ g exhibiting a clear nonlinear dependence as a function 
of the order q of the structure function. Rather than 
showing this standard looking multifractal spectrum, we 
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show in Fig. [3] another striking signature of multifrac- 
tality discovered first in empirical data j2^: the MRW 
predicts a continuous spectrum of exponents rj(s) for the 
relaxation of the volatility E[a 2 (t)\s] ~ t' 71 ^ from a lo- 
cal peak as a function of its amplitude cx e s given by 




o 1 = 8 
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realizations. The critical point of the Ising model is as- 
sociated with a critical value K c for the average coupling 
coefficient. Close to this value, agents organize spon- 
taneously within clusters of similar opinions, which be- 
come very susceptible to small external influences, such 
as a change of news. This may explain the occurrence of 
crashes as argued previously |23|. Intuitively, the critical 
slowing down well-known to characterize the proximity to 
the critical Ising point can explain the long-term memory 
of the volatility while the almost absence of correlation 
of the returns themselves is ensured by the impact of the 
news and the random idiosyncratic decisions. 

As a bonus, we have discovered that this simple model 
exhibits a rich multifractal structure, diagnosed not only 
by the standard convexity of the exponents of the struc- 
ture functions but also by a continuous spectrum of p ower 
law response functions to endogenous shocks [22l |24| . 
To our knowledge, this is the hrst nonlinear model in 
which such clear distinction is documented quantita- 
tively, based on a bottom-up self-organization. In con- 
trast, the multifractal random walk which has provided 
the theoretical predictions used here is a descriptive phe- 
nomenological mod el. 



FIG. 3: (Color online) For a given relative log-amplitude s of 
a local burst of volatility occurring at some time t a , we trans- 
late and superimpose all time series starting at those times 
of local bursts of the same amplitude s. Averaging over these 
time series of volatility obtains the average conditional relax- 
ation function of the volatility i?[a 2 (t)|s] ~ t~ v( - s ^ following a 
local burst of volatility of amplitude oc e s . The inset shows 
the average normalized conditional volatility _B[(j 2 |s]/_E[cr 2 ] as 
a function of the time after the local burst of volatility for dif- 
ferent log-amplitudes s. The figure shows the exponents rj(s) 
measured as the slopes of the curves in the inset for r = 1, 2, 4, 
and 8. 

We have outlined an Ising model with imitation be- 
tween agents, their influence by external news and the 
impact of their private information, which describes the 
"digestion" of external news by the collective behavior of 
the population of traders to create time series of returns 
presenting long-range memory in the volatility and mul- 
tifractal properties, in agreement with empirical data. 
This formulation has allowed us to test within the same 
model the hypothesis that agents are rational versus ir- 
rational. The empirical stylized facts of financial stock 
markets have been found only when agents misinterpret, 
or mis-attribute the source of the prediction of returns or 
are over-confident ({3 > 0). We can interpret our results 
by saying that, conditioned on their role of reflecting the 
stock market, the news serve as the substrate for fostering 
social interactions and reinforcing herding. Technically, 
the stylized facts in this regime result from the fact that 
the model operates around the critical point of the cor- 
responding Ising model, with coupling coefficients which 
are time-dependent and endowed with a memory of past 
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